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Introduction 

Desert surfaces are probably one of the most sta- 
ble of the Earth’s natural targets for remote sensing 
measurements. When viewed from space these sur- 
faces can, in effect, represent the solar diffuser plates 
or thermal blackbodies generally needed for in-flight 
instrument calibration. As an example, the bright 
gypsum sands and the nearby dark lava flows at the 
White Sands National Monument, New Mexico, are 
occasionally used as relative comparison targets for 
satellite sensors. St ay lor (1986) has discussed the se- 
lection and directional models of several deserts to be 
used as validation targets for three Earth Radiation 
Budget Experiment (ERBE) missions. These targets 
will serve as independent calibration sources in the 
event that the onboard calibration systems degrade 
or do not function as designed. 

The bidirectional reflectance properties of one of 
these deserts, the Saudi Arabian desert, was investi- 
gated by J. M. Davis (1982) during the Summer Mon- 
soon Experiment (Summer Monex) by using high- 
altitude aircraft at a time of large solar zenith an- 
gles. A comparison of the bidirectional reflectance 
determined from his data with that of Salomonson 
(1968) for low-altitude measurements of the White 
Sands desert showed significant differences between 
the two plots. (Davis reported an rms difference 
of 0.43.) The expectation was that the shape of the 
bidirectional reflectance curves should have been sim- 
ilar, even considering the differences in reflectances 
of deserts. Davis attributed these discrepancies to 
the forward scattering of the dust-laden atmosphere 
prevalent during Summer Monex. The three optical 
parameters related to the aerosol concentration and 
needed to compute the atmospheric effects are op- 
tical depth, single-scattering albedo, and scattering 
phase function. The relative importance of aerosol 
scattering and absorption in remote sensing has been 
discussed by Fraser and Kaufman (1983), and the ef- 
fects of surface reflectances on the radiation emerging 
from the top of the atmosphere have been presented 
by Fitch (1981). A study of the reflectance measure- 
ment errors produced by uncertainties in the input 
model parameters has been made by Bowker et al. 
(1983). 

The present paper is concerned with (1) a mod- 
eling of the influence of atmospheric aerosols on 
the remote sensing of bidirectional reflectances of 
desert targets, (2) a modeling of the effects of sur- 
face shadowing, and (3) a comparison of the results 
of these models with the flight results of Davis. This 
study was done only for large solar zenith angles be- 
cause the effects of aerosols and shadows are more 
pronounced at these angles and, more importantly, 


Davis’ data were obtained for solar zenith angles 
of 70° to 80°. Although it is possible to approxi- 
mate Davis’ results, it is felt that a surface reflectance 
model with a smaller specular component would have 
permitted using a more realistic set of atmospheric 
parameters in the simulations. 
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Volz turbidity parameter 

direct component of surface irradiance, 

W-m - 2 

sum of direct and diffuse components 
of surface irradiance, W-m -2 

cumulative frequency of daily average 
turbidity for one of three representa- 
tive distributions, percent 

forward/backward scattering ratio 

ratio of height of repetitive surface 
feature to wavelength 

altitude, km 

surface radiance resulting from diffuse 
irradiance component, W-m“ 2 -sr _1 

path radiance between surface and 
sensor, W-m _2 -sr _1 

surface radiance resulting from direct 
and diffuse irradiance components, 
W-m _2 -sr -1 

complex index of refraction 

real component of index of refraction 

imaginary component of index of 
refraction 

phase function 

atmospheric transmission from surface 
to sensor 

meteorological visibility, km 
visual range, km 
wavelength, fim 
scattering angle, deg 
sensor zenith angle, deg 
solar zenith angle, deg 
surface reflectance at nadir 
surface background reflectance 
optical depth 

azimuth angle (relative to Sun direc- 
tion), deg 


u 0 single-scattering albedo 

Abbreviations: 

AFGL Air Force Geophysics Laboratory 

BRF bidirectional reflectance factor 

ERBE Earth Radiation Budget Experiment 

ERIM Environmental Research Institute of 

Michigan 

Monex Monsoon Experiment 

Norm. refl. normalized reflectance 


rms root mean square 

Desert Reflectance Model 

The spectral reflectance of a collection of repre- 
sentative desert surfaces can be found in a report by 
Bowker et al. (1985). Most of these data were ob- 
tained with nadir looking, narrow-field-of-view sen- 
sors during small solar zenith angles. Ashburn and 
Weldon (1956) measured the reflectance of desert sur- 
faces by using a sensor with a 2tt steradian field of 
view; when these measurements are normalized by 
the total radiation incident on the surface, the results 
are referred to as albedo, an important term in clima- 
tological studies. Surface reflectances in the visible 
region of the spectrum varied from less than 0. 1 for 
dark basalt to greater than 0.7 for salt beds. In this 
investigation, then, reflectances of 0.1, 0.4, and 0.7 
will be taken to represent the average and typical 
range of values for desert environments. 

The reflections from most natural surfaces are 
anisotropic, particularly at large solar zenith an- 
gles. Kimes et al. (1985) have measured the bi- 
directional reflectance factors for several natural sur- 
faces of Northern Africa, ranging from plowed fields 
to salt plains. Salomonson (1968) determined the 
reflectance of a dry desert lake bed in California 
and the white gypsum sands near Alamogordo, New 
Mexico. The white sands were the most isotropic of 
the surfaces he studied. The normalized reflectances, 
i.e., the ratio of reflectance at the sensor zenith an- 
gle 0 r to the reflectance at nadir, of three desert sur- 
faces from Salomonson (1968) and Kimes et al. (1985) 
are shown in figure 1 for the principal plane. (The 
fundamental viewing geometry for this investigation 
is illustrated in figure 2, where the principal plane, 
defined as being normal to the surface and passing 
through the Sun and the origin of the reference sys- 
tem, is shown.) Note that the minimum reflectance 
of the white sands is located near nadir; these curves 
will be discussed in more detail later. The normalized 
bidirectional reflectance factor (BRF) of the white 
sands, for a solar zenith angle of 76° to 82°, is shown 


in figure 3; this is also referred to as an anisotropic 
reflectance diagram. Since Davis used Salomonson’s 
data in his paper, these data will be used here as 
an input in the radiative transfer calculations. For 
those interested in albedo, a technique for inferring 
the albedo of the Earth’s surface from multiple off- 
nadir view angles has been presented by Kimes and 
Sellers (1985). 

Two factors contributing to anisotropy in sur- 
face reflectance are surface roughness and shadow- 
ing. The importance of shadowing in bidirectional 
reflectances has recently been studied by Cooper 
and Smith (1985) and Norman et al. (1985). Both 
Salomonson (1968) and Kimes et al. (1985) at- 
tributed the minimum reflectance in the extreme 
forward-scattering direction of rough surfaces, such 
as soils and plowed fields, to shadowing. Insofar as 
the bidirectional data (fig. 3) selected for the present 
study appear to be relatively free of such effects, an 
independent model for shadowing is incorporated in 
the event that the surface includes such features as 
coarse regolith, dunes, and hills. Figure 4 illustrates 
the geometry of the surface roughness model in the 
principal plane where a terrain feature of height H 
has a periodic occurrence with horizontal wavelength 
A. The ratio H/X determines the percent of surface 
shadowing for a given solar zenith angle. The nomi- 
nal value of H/X — 0.134 used in this study gives a 
sunlight /shadow ratio of 1 at nadir for a solar zenith 
angle of 75°. For measurements out of the principal 
plane, shadowing is dependent on cos(ip). 

Radiative Transfer Program 

Radiative transfer determinations were made 
with two existing computer programs. A version 
of the ERIM radiative transfer model developed by 
Turner (1974) was used to account for single- and 
multiple-scattering effects in Rayleigh scattering and 
aerosol scattering as well as to account for ozone ab- 
sorption over a wavelength range of 0.27 to 2.20 /zm. 
The AFGL LOWTRAN 5 radiative transfer model, 
described by Kneizys et al. (1980), accounted for 
zero-order scattering and gaseous absorption effects 
over a wavelength range of 0.25 to 28.5 /zm. The 
spectral interval of interest for this investigation was 
limited to 0.55 to 0.85 /zm; all the data presented 
here are for the representative central wavelength 
of 0.7 /zm. 

In our version of the ERIM program, the target 
reflectance p is surrounded by a uniform background 
reflectance /?{,. The two reflectances are made equal 
in this study. The solar zenith angle 0 S was nomi- 
nally set at 75° to agree with Davis’ measurements, 
whereas the sensor nadir angle 0 r and the relative 
azimuth angle ip were varied at regular intervals. 


2 



Because of the assumed hemispheric symmetry of the 
bidirectional reflectance about the principal plane, 
the relative azimuth angle range simulated was lim- 
ited to 0°to 180°. The altitude h was set at 10 km 
to agree with the altitude of Davis’ measurements. 

The molecular effects (Rayleigh scattering) can 
be easily accounted for with a standard atmosphere, 
but the aerosol effects are variable and require the 
atmospheric conditions to be specified for each calcu- 
lation. Three optical parameters are needed for sim- 
ulation of the atmospheric effects: they are the opti- 
cal depth r, which determines the extinction of light; 
the single- scattering albedo which is that fraction 

of attenuation due to scattering; and the scattering 
phase function P(0 ), which gives the probability of 
scattering in a particular direction. 

Four values for the single-scattering albedo were 
investigated: 0.96, 0.81, 0.65, and 0.50. The 0.96 
value is representative of clear atmospheric condi- 
tions. Charlock and Sellers (1980) established a value 
of 0.81 for an atmosphere that was neither warm- 
ing nor cooling. The 0.65 value is a lower limit for 
dust-like material (Whitlock et al. 1985). The low- 
est value, 0.50, is more typical of an urban aerosol 
(Whitlock et al. 1985) or soot (Hansen et al. 1980). 

Figure 5 is a plot of the scattering phase functions 
for four different aerosols. The one labeled “IM.01L” 
was provided with the ERIM radiative transfer model 
and represents a typical land aerosol. It was used 
in all the calculations, whereas the others were se- 
lected for sensitivity studies only. (IM.01L signifies 
that the imaginary part of the refractive index (i.e., 
absorption coefficient) is equal to 0.01, and L signi- 
fies land aerosol. Model No. 25 is a model from the 
LOWTRAN 5 set by Kneizys et al. (1980).) An im- 
portant aspect of the curves is the forward/backward 
scattering ratio //£>; a high ratio implies a higher 
atmospheric scattering component of radiance when 
observation is in the direction of the Sun than when 
observation is in the direction opposite to the Sun. 

In the ERIM model, the effect of variability in 
the aerosol loading in the atmosphere is simulated 
indirectly by tying the aerosol optical thickness to 
the surface visual range VR through a built-in rou- 
tine. The routine scales the amount of aerosol below 
an altitude of 5 km to an estimate of the surface vi- 
sual range by decreasing the aerosol number density 
exponentially from the surface up to an altitude of 
5 km; above 5 km a constant background amount 
(independent of surface VR) is used. The common 
meteorological visibility V is related to visual range 
by VR = (1.3 ± 0.3)V (Kneizys et al. 1980). In our 
simulation VR was set equal to (1.3 -I- 0.3)V\ 


A more precise and meaningful way of simulat- 
ing aerosol scattering is to input such optical depth 
values directly rather than to rely on> the approx- 
imate visibility relationship. Flowers et al. (1969) 
have established the cumulative frequencies of daily 
average turbidities for three representative distribu- 
tions (rural, suburban, and urban) based on the Volz 
aerosol turbidity parameter B. The optical thickness 
r is related to turbidity by the relation e~ T = 10“ 5 . 
Thus, r = 2.30259H. The aerosol turbidity param- 
eter B is determined by looking at the Sun with a 
standard photometer, at a wavelength of 0.5 //m, 
and correcting for the optical thickness contributions 
due to molecular scattering and ozone absorption. 
Table I gives the relations between turbidity, visual 
range, and visibility for the three distributions, for 
the 5th, 50th, and 95th percentiles. (In explanation, 
it should be noted that the 50th percentile is equiv- 
alent to the average condition. The 5th percentile 
corresponds to a quite clear atmosphere, exceeded in 
turbidity 95 percent of the time; the 95th percentile 
corresponds to a very hazy atmosphere, exceeded in 
turbidity only 5 percent of the time.) A rural at- 
mosphere was used in this study, since it is believed 
most typical of desert environments; thus, the tur- 
bidity varied from 0.035 to 0.125 (at 0.5 //m) in the 
radiative transfer calculations. Turbidity values at 
other wavelengths are extrapolated. 

Table I. Values of Turbidity B, Visual Range VR, 
and Visibility V for Various Percentile Values, 

From Flowers et al. for Three Representative 
Cumulative Frequency Distributions 


[Estimated VR = (1.6)V] 


Parameter 

Rural 

Suburban 

Urban 

5th percentile 




Turbidity 

0.035 

0.038 

0.055 

Visual range, km 

100 

97 

73 

Visibility, km 

76 

74.6 

56.2 

Visibility, mi 

47.8 

46.3 

34.9 

50th percentile 




Turbidity 

0.063 

0.094 

0.175 

Visual range, km 

56 

37 

14.5 

Visibility, km 

43.1 

28.5 

11.2 

Visibility, mi 

26.8 

17.7 

6.9 

95th percentile 




Turbidity 

0.125 

0.275 

0.450 

Visual range, km 

26 

7.6 

4 

Visibility, km 

20 

5.8 

3.1 

Visibility, mi 

12.4 

3.6 

1.9 


The radiative transfer computer program gives 
the radiance at the sensor for each set of input values. 
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However, when part of the surface is in shadow a 
separate calculation must be made. The radiance 
at the sensor from the shaded area is the result of 
scattered sunlight and is equal to 

L(indir) = ([E(tot) - 25(dir)] pT / 7r) + L(path) 

Therefore, the total radiance at the sensor is a com- 
bination of L(tot) and L(indir), according to the sun- 
light/shadow ratio. 

Results 

Radiance determinations were made for an alti- 
tude of 10 km and a solar zenith angle of 75° for each 
of 36 possible combinations of surface reflectance, 
single-scattering albedo, and turbidity percentile val- 
ues previously discussed. The IM.01L aerosol model 
was used. The sensor zenith angle and the relative 
azimuth angle were varied in 15° and 22.5° incre- 
ments, respectively, which resulted in 46 calculations 
per set of input conditions. The radiance values were 
normalized to nadir radiance, plotted on a hemi- 
spheric grid displaying 9 r and tp, and contoured at 
regular intervals. A representative set of these nor- 
malized bidirectional reflectance plots is shown in fig- 
ures 6 through 11 for p = 0.4 and 0 9 = 75°. In 
figures 6 through 9, F(B) = 50 percent (i.e., the 
50th percentile or average turbidity condition) and 
oj 0 varies from 0.96 to 0.50. For u 0 ~ 0.81, two addi- 
tional plots are shown for F(B) values of 95 percent 
and 5 percent in figures 10 and 11. 

To compare the off-nadir reflectance variations re- 
sulting from changes in the single-scattering albedo 
and turbidity values, normalized reflectances in the 
principal plane have been plotted for surface re- 
flectances of 0.1, 0.4, and 0.7 in figures 12, 13, and 14, 
respectively. The upper part of each figure is for 
u 0 = 0.96, and the lower part is for c o Q = 0.50. Each 
figure displays five curves: the two extreme condi- 
tions ( F(B ) = 95 percent and 5 percent) without sur- 
face shadowing and with shadowing (H / A = 0.134), 
and the surface reflectance from Salomonson (1968) 
that was used in the computations. 

Figure 15 was redrawn from Davis (1982) and 
shows contours of differences between the normalized 
reflected radiances from his study and the same 
quantity from the data of Salomonson (1968) for a 
solar zenith angle of 70° to 80°. Positive values 
indicate that brighter features were measured by 
Davis. Difference values (referred to as A norm, refl.) 
were extracted for the principal plane and used as a 
basis for comparison in the next four figures, which 
show the differences between calculated values and 
the Salomonson data for various input conditions. 


Figure 16 shows the effect on A norm. refl. of 
variations in u 0 and p for a solar zenith angle of 75°, 
an F(B) of 5 percent, and shadowing. In figure 17 
the solar zenith angle is varied from 70° to 80° with 
constant input conditions. Figures 18 and 19 are 
for a solar zenith angle of 70°, a single-scattering 
albedo of 0.50, and aerosol model No. 25, which has 
a low forward/backward scattering ratio. Figure 18 
illustrates the influence of shadowing variations for 
F(B) values of 5 percent and 50 percent, for p = 0.4, 
while figure 19 gives the same for p = 0.7. 

Discussion 

Aerosols 

Figures 6 through 9 illustrate the effect of the 
single-scattering albedo on reflectance measure- 
ments at high altitude. With no intervening atmos- 
phere, these hemispheric plots should look like the 
Salomonson data in figure 3, which was used as the 
input to the program. Scattering and absorption, 
however, have caused a distortion in the reflectance, 
particularly along the principal plane. The highest 
scattering value (cj q = 0.96) introduced the greatest 
error in the measurement because of the increased 
path radiance at high zenith angles. In figures 10 
and 11, where the single-scattering albedo is con- 
stant and the turbidity is varied, the higher turbidity 
produced the most distortion in the reflectance pat- 
tern because of the increased path radiance. There 
are other factors that could also influence these re- 
flectance plots, but they will be discussed later when 
comparisons are made with the data from Davis. 

Notice that the reflectance plots, although 
distorted by the atmosphere, still resemble the 
Salomonson data. The reflectance data from Davis 
presented in figure 15 shows more detail, especially 
at nadir, than the original Salomonson data. These 
features may be real or an artifact of the measure- 
ments. Since Davis collected data in only 12 off-nadir 
directions at angles of either 30° or 60°, and then per- 
formed a harmonic analysis to generate the final bi- 
directional reflectance plot, the variations at nadir 
are considered to be artificial. On the other hand, 
they may have been suppressed in the data presented 
here by the scattering phase function entered into the 
program and/or the simplifications used in the radia- 
tive transfer calculations within the program. 

Because of the subtle differences in the normal- 
ized reflectances when the atmospheric and surface 
conditions are varied, comparisons of the data are 
best shown in the principal plane. From the plots 
displayed in figures 12 through 14, four general ob- 
servations can be made. First, an increase in tur- 
bidity increases the normalized reflectances at high 
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zenith angles, particularly when looking in the direc- 
tion of the Sun (ip = 0°). This increase is primarily 
due to the increased atmospheric path length relative 
to nadir and to the fact that the aerosol phase func- 
tion is peaked in the forward direction and to a lesser 
extent in the backward direction (see fig. 5), which 
causes the path radiance to be more pronounced. 
Second, a decrease in the single-scattering albedo 
tends to reduce the relative atmospheric effect at 
high zenith angles, which results in a lower normal- 
ized reflectance. This decrease is brought about by 
an increase in the aerosol absorption, which has a 
more pronounced effect over the longer path length. 
Third, when the surface reflectance is increased, the 
normalized reflectances at high zenith angles are re- 
duced. In this case the normalizing effect of the 
higher nadir radiance, which increases directly with 
surface reflectance, is responsible. And fourth, the 
surface shadowing effect tends to decrease the nor- 
malized reflectance in the ip = 0° direction and in- 
crease it in the ip = 180° direction. Since shadow- 
ing has reduced the radiances continuously from the 
ip = 180°direction to the ip = 0° direction, the net 
effect is to decrease the radiance relative to nadir in 
the ip = 0° direction and to increase the radiance 
relative to nadir in the ip = 180° direction. 

Although the above trends in the data were gen- 
erally anticipated, it is the magnitudes of the effects 
that are most important. Several of the normalized 
reflectance curves show a decrease in reflectance rela- 
tive to the surface in the ip — 0° direction. Since the 
Davis data also behave in this manner, it is meaning- 
ful to note those conditions that enhance this effect. 
It is seen that low turbidity, low single-scattering 
albedo, high surface reflectance, and shadowing are 
generally favorable to a reduction in the normalized 
reflectance. These variables, along with several oth- 
ers, will next be manipulated in the input parameters 
in an attempt to reproduce the results of Davis. 

Davis Flight Data 

In figures 16 through 19 several input param- 
eters were varied and the results compared with 
those of Davis. Figure 16, which illustrates the ef- 
fect of varying surface reflectance for two single- 
scattering albedo conditions, shows that the more 
absorbing atmosphere and the higher reflectance sur- 
face produced the most favorable results. In figure 17 
the most favorable turbidity (5 percent) and single- 
scattering albedo (0.50) conditions were held con- 
stant, and the solar zenith angle was varied. The 
zenith angle of 70° gave the best results in the ip = 0° 
direction but least favorable in the ip = 180° direc- 
tion because of the highly peaked scattering phase 
function. Decreasing the sensor-surface-Sun angle, 


by elevating the Sun, reduced the path radiance rel- 
ative to nadir, and thus depressed both ends of the 
curve. 

In figure 18 the effect of varying two new param- 
eters is illustrated. The aerosol phase function has 
been changed with the use of LOWTRAN 5 model 
No. 25 (which has a low f/b ) shown in figure 5, 
and the shadowing has been further increased to 
H/X = 0.2. With the solar zenith angle equal to 70°, 
the three shadowing values of 0.1, 0.134, and 0.2 give 
28, 37, and 55 percent shadow at nadir, respectively. 
Increasing the shadowing reduced the reflectance in 
the ip = 0° direction and increased it in the ip = 180° 
direction. The curve for H/X = 0.134 from figure 17 
lies between the two extremes in figure 18. Although 
the new aerosol model has changed the shape of these 
curves somewhat, the major influence has been to re- 
duce the extreme values in the ip = 0° direction. In- 
creasing the surface reflectance in figure 19 had only 
a small effect on the results. 

The variation in Davis’ data at nadir can be ap- 
proximated by varying the surface reflectance across 
the field of view: for instance, by letting p — 0.4 in 
the ip = 0° direction and p = 0.7 in the ip = 180° 
direction, with the nadir reflectance being an aver- 
age of the two values (i.e., 0.55). Although such a 
change in surface reflectance is physically possible, it 
does not seem reasonable for data collected during a 
solar zenith angle variation of 10°. As mentioned be- 
fore, this feature of Davis’ data is thought to be an 
artifact resulting from the mathematical harmonic 
analysis of limited data. 

It is clear in the simulation of Davis’ data that 
the values of the atmospheric parameters are not 
what one would expect during the summer monsoon 
season. Indeed, some of them are quite unrealistic, 
such as a very clear atmospheric condition with a 
low value of single-scattering albedo. In the following 
section, a discussion will be given concerning a viable 
range of environmental parameters for desert areas. 

Desert Environments 

There are several good sources of information 
concerning desert environmental conditions. How- 
ever, none of them give all of the parameters nec- 
essary to adequately describe the radiative transfer 
conditions. For the purpose of this discussion, six im- 
portant variables will be reviewed; they are aerosol 
distribution, optical depth, single-scattering albedo, 
scattering phase function, surface reflectance, and 
shadowing. 

Patterson and Grams (1981) investigated the 
haze over the Saudi Arabian peninsula and the 
Arabian Sea during the summer monsoon season and 
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found that the aerosol layer extended from the sur- 
face to an altitude of approximately 6 km over the 
desert regions. (It should be noted here that the 
Patterson and Grams data and the Davis data were 
collected 2 years apart.) Over both areas the highest 
concentrations were near the surface, and there was 
an indication of layering. There was not one well- 
mixed layer extending from the surface to the top of 
the haze, but multiple layers. In both of these cases, 
the concentration actually decreased at the lowest 
level measured, which implies that the source of the 
aerosol was not directly below the measurement lo- 
cation and suggests that transport from some distant 
source was important. The layering also suggests 
the possibility of multiple sources for the haze layers. 
Schutz and Sebert (1984) also found this transport 
condition to be characteristic of the Saharan dust 
layer. 

As previously mentioned, the Turner program 
used in this study assumes a well-mixed aerosol layer 
with an exponential decrease from the surface to 
an altitude of 5 km. This aspect of the program 
is considered to have a minor influence on the re- 
sults, however, since the aerosol effects are usually 
variable with time and distance and are difficult to 
model. A more important parameter is the optical 
depth, which integrates the effect of the total aerosol 
burden. 

In a review of global atmospheric aerosols, Toon 
and Pollack (1976) referenced an optical depth over 
Israel of r = 0.24 when dust is blown from the North 
African deserts. For the Saudi Arabian peninsula, 
Patterson and Grams (1981) described the visibility 
as less than 10 km and assumed an optical depth of 
approximately 1.5. Although these two conditions at 
first glance do not seem to be compatible with table I 
(that is, a visibility of 3 km represents an optical 
depth of only 1.2), it should be remembered that the 
aerosol density actually decreased near the surface, 
which permitted a higher visibility than suggested by 
the optical depth. Endemann et al. (1984) developed 
an aerosol model that was thought to be typical of a 
heavy haze in desert areas of the world, especially in 
spring and summer. The continental-type boundary 
layer extended up to 6 km and had an optical depth 
of 3.3. This value exceeds the boundaries of the 
turbidity models that were established by Flowers 
et al. (1969) for the United States and were used in 
this study. 

The optical effects of the aerosols, which include 
the single-scattering albedo and the scattering phase 
function, are determined by the complex index of re- 
fraction, the size distribution, and the shape of the 
aerosol particles. The complex index of refraction 
n = n ve — in im consists of a real part n re , which is 


the ratio of the velocity of light in a vacuum to the 
velocity of light in the material, and an imaginary 
part rii m , which is an absorption parameter charac- 
terizing the material. Patterson et al. (1977) and 
Patterson (1981) have evaluated the optical proper- 
ties of crustal aerosols and found for visible wave- 
lengths that n re lies between 1.5 and 1.6 and n im lies 
between 0.001 and 0.1. The real part controls scat- 
tering whereas the imaginary part determines aerosol 
absorption. By assuming that the particles are spher- 
ical, the Mie theory (see Bohren and Huffman, 1983) 
can be used to determine the interaction of the par- 
ticles with the radiation field. A number of authors 
have constructed aerosol models and used the Mie 
theory to define the optical effects. 

Braslau and Dave (1973a and b) have studied 
aerosols with an imaginary refractive index of 0.01 
and found that uj 0 was about 0.9. Patterson and 
Grams (1981) considered the Saudi Arabian desert 
haze with an optical depth of 1.5 and determined 
an atmospheric absorption of about 15 percent; i.e., 
(jj 0 — 0.85. Minnis and Cox (1978) found that 
dust in the Saharan air layer could cause increases 
in the fractional absorption of the atmosphere from 
4 percent to 19 percent. 

In their study of the Arabian peninsula aerosol 
layer, Patterson and Grams noticed a highly ab- 
sorptive background material, assumed to be carbon, 
above the haze layer. The effect this material has on 
the optical properties of the aerosol layer will depend 
on its concentration. The desert aerosol model of En- 
demann et al. (1984) contained 1 percent soot and 
had an u) 0 of 0.88. Their urban-industrial haze model 
was composed of 22 percent soot with u Q = 0.64. 
Bohren and Huffman (1983) calculated the effect of 
an addition of 1 percent soot to an aerosol model with 
two different techniques. In the first case the small 
soot particles were mixed with the larger (1.5 /mi) 
particles, and in the second case the soot was imbed- 
ded in the larger particles. The imbedding technique 
enhanced the effectiveness of the absorption of the 
soot and decreased w 0 from about 0.9 to 0.75. This 
method of mixing could be very important in the 
event that volatile aerosols are present and not re- 
covered in the sampling procedure. 

The aerosol-scattering phase function, which is 
wavelength dependent, usually has a maximum in 
the forward direction, a lesser maximum in the back- 
ward direction, and a minimum between 90° to 120°. 
According to Barteneva (1960) the phase function is 
nearly flat for a very clear atmosphere and gradu- 
ally becomes more forward scattering as the turbid- 
ity increases. Kneizys et al. (1980) have developed 70 
aerosol-scattering phase functions that are consid- 
ered to be representative of various aerosol models 
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and wavelengths. Model No. 25, used in this re- 
port, has a slightly higher forward/backward scatter- 
ing ratio than a very clear atmosphere but much less 
than the phase functions used in their tropospheric 
models. (It is important to note that when multi- 
ple scattering becomes important, such as in a dense 
cloud, the surface tends to become more diffuse, 
which reduces the effect of a high / /b ratio.) In the 
visible region, the forward/backward scattering ratio 
of common land aerosols can vary from about 40 to 
over 400. 

The reflectance of the White Sands target was de- 
scribed as the most isotropic target investigated by 
Salomonson. For higher solar elevation angles, the 
principal plane reflectance showed only a slight in- 
crease in the ip = 180° direction and essentially no 
increase in the ip = 0° direction. As the solar zenith 
angle increased, the reflectance increased somewhat 
in the ip = 180° direction and even more in the 
ip = 0° direction. The increase in the ip = 0° direc- 
tion is taken to be due to specular reflectance. For 
those desert surfaces with vegetation, the reflectance 
increased back toward the Sun (in the ip = 180° di- 
rection) and decreased away from the Sun. This oc- 
currence can be described qualitatively as a shadow 
effect that appears to be related to the size of the 
elements creating the reflection. In the reflectances 
of the dry desert lake bed and the salt plain sur- 
faces, shown in figure 1, the small humps between 
nadir and the solar zenith angle are related to the 
“hot spot”, i.e., the point where the Sun-sensor line 
intersects the surface. In laboratory measurements 
the hot spot was located at the solar zenith angle, 
whereas the displacement in the field was thought 
to result from a more diffuse illumination condition 
(Salomonson, 1968). 

Although shadowing is an integral part of the 
reflectance mechanism of most natural surfaces, it 
is essential to consider the influence of larger fea- 
tures not generally accounted for in field measure- 
ments. Rocks, boulders, dunes, hills, and even moun- 
tains can be important to the scene reflectance when 
viewed from high altitude. Landsat images of the 
Arabian peninsula (Short et al., 1976) reveal a va- 
riety of dune types and mountainous terrain where 
shadowing is clearly evident even with a relatively 
small solar zenith angle. In a review of regional land- 
forms, Short and Blair (1986) noted eolian features 
with heights of 50 to 150 m and wavelengths of 500 m 
to 3 km; these could account for about one-half of the 
shadowing in the nominal case (H / A = 0.134) used 
in this study. 

In closing the discussion on desert environments, 
it is apparent that the uncertainties in the scatter- 
ing function, single-scattering albedo, and turbidity, 


even under favorable conditions, preclude the selec- 
tion of a most probable set of these parameters for 
the radiation transfer study. However, one parame- 
ter that has a significant effect on Davis’ data and 
is independent of these atmospheric parameters is 
the surface reflectance. The reflectance model cho- 
sen contains a specular reflectance component that 
has constrained the choice of atmospheric turbidity 
for the simulation of Davis’ flight data. Salomon- 
son remarks that the most likely principal plane re- 
flectance curve for sand would have a minimum near 
nadir and would increase toward the Sun. Remov- 
ing the specular component from Salomonson’s data 
and introducing a shadowing factor would permit a 
more reasonable choice of aerosol conditions for the 
comparison with Davis’ results. 

Concluding Remarks 

A review of atmospheric conditions for desert ar- 
eas has revealed that environmental parameters can 
be highly variable, which makes it difficult to estab- 
lish an average condition for radiative transfer sim- 
ulations. In this study concerning the influence of 
aerosols on the remote sensing of surface reflectances, 
a rural atmospheric condition (Flowers et al., 1969) 
was assumed to be typical of North American deserts. 
A desert reflectance for a high solar zenith angle, 
taken from Salomonson (1968), was used as the sur- 
face reflectance model. With a suitable choice of in- 
put values for the modeling program, it was possi- 
ble to approximate the flight results of Davis (1982) 
for the Saudi Arabian desert. The assumed low tur- 
bidity conditions did not seem reasonable, however, 
particularly for that region of the globe during the 
summer monsoon season. It was concluded that a dif- 
ferent surface reflectance model, e.g., one without a 
specular component and with some shadowing, would 
have permitted a more realistic choice of modeled at- 
mospheric conditions and should have given a closer 
match to the observed reflectance curves. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
August 5, 1987 
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Figure 1. Normalized reflectances in the principal plane for three typical desert surfaces: white sand and dry 
desert lake bed (Salomonson, 1968), and salt plain (Kimes et al., 1985). The solar zenith angle is shown 
for each plot. 
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Figure 4. Illustration of the shadowing model used in the radiative transfer calculations for the principal plane; 
the height H of a repetitive surface feature with wavelength A is shown. The feature is assumed to be 
continuous and perpendicular to the paper. 
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Figure 5. The normalized scattering phase functions for four aerosol models used in the radiative transfer 
calculations. 
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Figure 10. Normalized BRF plot for desert surface with p = 0.4, oj q = 0.81, F(B) = 95 percent, and 0 9 = 75°. 
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Figure 11. Normalized BRF plot for desert surface with p = 0.4, oj 0 = 0.81, F(B) = 5 percent, and 6 S = 75° 
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Figure 12. Normalized reflectances in the principal plane for p = 0.1, 0 S — 75°, and u> 0 = 0.96 and 0.50. The 
surface reflectance is shown for comparison in figures 12 to 14. 



Norm. 

refl. 




Figure 13. Normalized reflectances in the principal plane for p — 0.4, 9 S — 75°, and oj 0 — 0.96 and 0.50. 
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Figure 14. Normalized reflectances in the principal plane for p = 0.7, 9 S = 75°, and lo q — 0.96 and 0.50 
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Figure 15. Plot of differences between the reflectance model of Salomonson in figure 2 and the flight results of 
Davis (1982). (Redrawn from Davis, 1982.) Positive values indicate that brighter features were measured 
by Davis. 
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